Phase diagrams of colloidal spheres with a constant zeta-potential 
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We study suspensions of colloidal spheres with a constant zeta-potential within Poisson-Boltzmann 
theory, quantifying the discharging of the spheres with increasing colloid density and decreasing salt 
concentration. We use the calculated renormalized charge of the colloids to determine their pair- 
wise effective screened-Coulomb repulsions. Bulk phase diagrams in the colloid concentration-salt 
concentration representation follow, for various zeta-potentials, by a mapping onto published fits of 
phase boundaries of point- Yukawa systems. Although the resulting phase diagrams do feature face- 
centered cubic (fee) and body-centered cubic (bee) phases, they are dominated by the (re-entrant) 
fluid phase due to the colloidal discharging with increasing colloid concentration and decreasing salt 
concentration. 



I. INTRODUCTION 

Charged colloidal particles suspended in a liquid elec- 
trolyte are interesting soft-matter systems that have gen- 
erated fundamental as well as industrial attention for 
decades. pQ Understanding the stability and phase be- 
haviour of these systems as a function of colloid con- 
centration and ionic strength is an important theme in 
many of these studies. A key role is played by the 
electrostatic repulsions between the colloidal spheres, 
which are not only capable of stabilising suspensions 
against irreversible aggregation due to attractive Van 
der Waals forces, [2] but are also the driving force for 
crystallisation, [3] provided the surface charge on the col- 
loids is high enough and the range of the repulsions long 
enough. [TH3] The classic theory that describes the elec- 
trostatic repulsions between charged colloidal particles in 
suspension goes back to the 1940's, when Derjaguin, Lan- 
dau, Verwey, and Overbeek (DLVO) found, within lin- 
ear screening theory, that suspended spheres repel each 
other by screened-coulomb (Yukawa) interactions. |H [5] 
The strength of these repulsions increases with the square 
of the colloidal charge, and they decay exponentially with 
particle-particle separation on the length scale of the 
Debye screening length of the solvent. [6] This pairwise 
Yukawa form is a corner stone of colloid science, and can 
explain a large number of observations. P~H3] For instance, 
the experimentally observed crystallisation of charged 
colloidal spheres into body-centered cubic (bec) and face- 
centered cubic (fee) phases upon increasing the colloidal 
packing fraction at low and high salt concentrations, [TJ- 
110] respectively, is in fair agreement with simulations 
of Yukawa systems. |lllfl4] Interestingly, in these simu- 
lation studies, as well as in many other studies. [T5HT8] 
the charge of the colloids is assumed to be independent 
of the colloid density and the salt concentration. 

The constant-charge assumption was argued to break 
down, however, in some recent studies where the electro- 
static repulsions were argued to be reduced with increas- 
ing colloid concentration. Biesheuvel. [19] for instance, 
argues that experimental equilibrium sedimentation- 



diffusion profiles of charged silica spheres in ethanol at 
extremely low salt concentrations [20] are better fitted 
by a charge-regulation model than by a constant-charge 
model. [H] More recent evidence for a concentration- 
dependent colloidal charge stems from re-entrant melting 
and re-entrant freezing observations of PMMA spheres 
in a solvent mixture of cis-decaline and cyclohexyl bro- 
mide, i.e. the phase sequence upon increasing the col- 
loid concentration is fluid-crystal- fluid-crystal. [2"2"j In ad- 
dition, direct force measurements between a single pair 
of colloidal PMMA spheres in hexadecane, a pair that is 
part of a triplet, and a pair that is part of a multiplet have 
very recently revealed a significant reduction of the force 
with increasing number of neighbouring particles. |23j In- 
terestingly, in the three experiments of Refs.[3Ul [S3] 
the solvent is a nonpolar medium. 



In fact, the experimental findings of Ref . [23] could 
well be interpreted and explained in terms of constant- 
potential boundary conditions on the colloidal surfaces, 
rather than the more usual constant-charge assump- 
tion. The present article addresses the consequence of 
constant-potential boundary conditions for the packing 
fraction-salt concentration phase diagram of Yukawa sys- 
tems by calculating the colloidal charge and the effective 
screening length for various zeta-potentials as a func- 
tion of salt- and colloid concentration. Perhaps surpris- 
ingly, such a study has not yet been performed. In the 
case of high zeta-potential this requires nonlinear screen- 
ing theory, and hence the renormalized rather than the 
bare colloidal charge determines the effective screened- 
Coulomb repulsions between the colloids. [13 H3 H3H5U] 
For this reason we use the renormalized charge through- 
out. We also compare our constant-potential calculations 
with those of an explicit charge-regulation model, [3D - f3l)] 
and conclude that their results are qualitatively similar, 
and even quantitatively if they are considered as a func- 
tion of the effective screening length. 
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II. MODEL AND THEORY 

We consider N colloidal spheres of radius a in a sol- 
vent of volume V, temperature T, dielectric constant e 
and Bjerrum length Xb = e 2 /eksT. Here e is the el- 
ementary charge and fcs the Boltzmann constant. The 
colloidal density is denoted by n = N/V and the pack- 
ing fraction by 77 — (47r/3)na 3 . The suspension is pre- 
sumed to be in osmotic contact with a 1:1 electrolyte of 
Debye length k _1 and total salt concentration 2p s . We 
are interested in suspensions of charged colloids of which 
the surface (zeta) potential ipQ rather than the charge 
Ze is fixed. We will show that this constant-potential 
condition mimics charge-regulation on the colloidal sur- 
faces fairly accurately. The first goal of this article is 
to calculate Z as a function of r\ for fixed dimensionless 
combinations na, a/ Xb, and 4>q = eipa/ksT. This result 
will then be used to quantify the effective Yukawa inter- 
actions between pairs of colloids, and hence the phase 
boundaries between fluid, face-centered cubic (fee) and 
body-centered cubic (bec) crystalline phases. 

In the actual suspension of constant-potential colloidal 
spheres, the charge distribution of each of the N col- 
loids will be distributed heterogeneously over its surface 
due to the proximity of other colloids in some directions. 
This leads to a tremendously complex many-body prob- 
lem that we simplify here by assuming a spherically sym- 
metric environment for each colloid, which is nevertheless 
expected to describe the average electrostatic properties 
realistically. Below we will calculate the electrostatic po- 
tential ip(r) at a radial distance r from a charged colloidal 
sphere at a given zeta-potential ipo, le. at a given value 
ip(a) = x/jq. The colloidal charge Ze then follows from 
Gauss' law 



Ze 

9 ' 



(i) 



where a prime denotes a radial derivative. 

We first consider a single colloid in the center of a 
Wigner-Seitz cell of radius R, such that the cell vol- 
ume equals the volume per particle, (47r/3)i? 3 = V/N, 
which implies R = a?/" 1 / 3 . The radial coordinate of 
the cell is called r. We write the ionic density pro- 
files for r € (a, R) as Boltzmann distributions p± (r) = 
p s exp(=F0(r)), with <p(r) = eip(r)/k B T the dimension- 
less electrostatic potential. Together with the Poisson 
equation V 2 = — 47tAb(/9 + (7-) — p_(r)), this gives rise to 
the radially symmetric PB-equation and boundary con- 
ditions (BC's) 

4>"{r) + VO) = K 2 sinh0(r), re(a,R); (2) 
r 

<P(a) = 0o; (3) 

<t>'{R) = o, (4) 

where a prime denotes a derivative with respect to r. 
Note that BC Q implies charge neutrality of the cell. 
Once the solution 4>{r) is found for given 77, Ka, and <fio, 



e.g. numerically on a radial grid, the colloidal charge Z 
follows from Eq. 0, which we rewrite in dimensionless 
form as 



ZX r 



= —a<f> (a). 



(5) 



From the numerical solutions that we will present below 
it turns out that Z decreases monotonically from a finite 
asymptotic \ow-r] (large-i?) value Zq to essentially at 
77 ~ 1 (or R ~ a). 

Within linear-screening theory at low packing fraction, 
where sinh <f> ~ <fi, the potential distribution can be solved 
for analytically, yielding <j){r) = 0oacxp[— n(r — a)]/r, 
such that the colloidal charge takes the asymptotic low-77 
and low-</>o value 



Z n X 



0*B 



(1 + K,a)<f>o- 



(6) 



In the appendix we show that the discharging effect with 
increasing 77, as found from the nonlinear screening the- 
ory discussed above, can also be approximated within 
linear screening theory, yielding 



Z(tj, na) — 



(ko) 



1 + 77/77* 



3(1 + ko)' 



(7) 



where 77* is a crossover packing fraction at which the 
colloidal charge has decayed to half its dilute-limit value 
Zo given in Eq.([6|. For typical numbers of experimental 
interest, e.g. a/Xs = 100 and na = 0.25, we then find 
Z = 12500 and 77* = 0.017. With O ~ 1 - 2, which 
corresponds to a surface potential of 25-50mV, we should 
expect a few hundred charges in the dilute limit and a 
significant charge reduction for 77 > 10 -2 . 

The constant-potential boundary condition that we 
employ here is supposed to mimic charge-regulation on 
the colloidal surface through an association-dissociation 
equilibrium of chargeable groups on the surface. Here 
we consider, as a typical example, the reaction SA^ 
S + +A~, where a neutral surface group SA dissociates 
into a positively charged surface group S + and a re- 
leased anion A - . The chemistry of such a reaction can 
be characterised by a reaction constant K such that 
[S + ][A~]/[SA]=iT, where the square brackets indicate 
concentrations in the vicinity of the surface where the 
reaction takes place. If we now realise that Z cx [S + ], 
we find for the usual case where [S + ] -C[SA] that Z oc 
1/[A~]. For the case that the released anion is of the 
same species as the anion in the reservoir, such that 
[A~] = p s exp[<fi(a)], we thus have 



Z = z exp(— <j>(a)), 



(8) 



where the prefactor z, which is a measure for the sur- 
face chargeability, [37] accounts for the chemistry, the 
surface-site areal density, and the total area of the surface 
between the colloidal particle and the electrolyte solu- 
tion. Note that Eq.([8]) relates the (yet unknown) colloidal 
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charge Z to the (yet unknown) zeta-potential 0(a), for 
a given z. A closed set of equations for charge-regulated 
colloids is obtained by combining the PB equation ^ 
with BC Q at the boundary of a Wigner-Seitz cell of 
radius R, with BC ^ replaced by 



a(f> (a) = — cxp(— (j){a)), 

a 



(9) 



for some given chargeability z. The resulting solution 
</>(r) gives the zeta-potential (j)(a) as well as the colloidal 
charge Z using Eq.(|8]). When comparing the constant- 
potential model with the ionic association-dissociation 
model, we will tune the chargeability z such that the 
low-77 results for Z coincide for both models. 

It is well known that nonlinear screening effects, in 
particular counterion condensation in the vicinity of a 
highly charged colloidal surface, reduce the effective col- 
loidal charge that dictates the screened-Coulomb inter- 
actions between the colloids. [HI OH I24H26] 138] The so- 
called renormalized colloidal charge, i?*,[3H] can be cal- 
culated from the electrostatic potential <fi(r) as obtained 
from the nonlinear PB equation by matching the numer- 
ically obtained solution at the edge of the cell to the ana- 
lytically known solution of a suitably linearized problem. 
By extrapolating the solution of the linearized problem 
to the colloidal surface at r = a, one obtains the effective 
charge by evaluating the derivative at r — a using Eq. ^ . 
Following Trizac et qZ. . [40] the renormalized charge Z* 
can be written as 



Z*X, 



tanh < 



((R 2 aR - 1) sinh[K(i? - a)} + 



k(R — a) cosh[K(i? — a)]) , 



(10) 



where the 'Donnan' potential is defined as 4>d = (f>{R)i 
i.e. the numerically found potential at the boundary of 
the cell, and where the effective inverse screening length 
is 



K = K\/ COsh (/)£,. 



(11) 



Note that Z* and R can be calculated for the constant- 
potential as well as the association-dissociation model in 
a Wigner-Seitz cell. 



III. EFFECTIVE CHARGE AND SCREENING 
LENGTH 

For both the constant surface potential (CSP) and 
the association-dissociation (AD) model discussed above 
we calculated the bare colloidal charge Z, the effec- 
tive (renormalized) charge Z* , and the effective inverse 
screening length R in the geometry of Wigner-Seitz cells. 
In Fig[l] we show Z\s/a (full curves) and Z*Xb/ci 
(dashed curves) as a function of packing fraction 77, for 
two screening constants for both the CSP model (black 
curves) and the AD model (red curves), in (a) for fixed 
zeta-potential cf>o — 1 and in (b) for 4>q = 5. In all cases 
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Figure 1: The bare colloidal charge Z (continuous black 
curves) and the renormalized charge Z* (dashed black 
curves), both in units of a/As (see text), as a function of 
the colloidal packing fraction rj for several screening param- 
eters na, for constant surface potentials (a) cj>o = 1 and (b) 
4>o =5. The red curves denote Z and Z* as obtained from the 
association-dissociation model, with the chargeability z cho- 
sen such that the surface potential in the dilute limit 77 — >■ 
equals 0o- 



the chargeability parameter z of the AD model is cho- 
sen such as to agree with the CSP model in the low- 
density limit 77 — > 0. The semi-quantitative agreement 
between the red and black curves for equal na is indica- 
tive of the reasonable description of charge-regulation by 
constant-potential BC's, certainly at low 77. At higher 
?7 the charges predicted by the AD model exceed those 
of the CSP model somewhat, which should not come 
as a surprise since the former interpolates between the 
constant-charge and the constant-potential model. The 
close agreement between Z and Z* for all Ka at (f>o = 1 
in Fig[lja) is also to be expected, since <f>Q — 1 is not 
far into the nonlinear regime. By contrast, deep in the 
nonlinear regime of 4>$ — 5, as shown in FigJIJb), there 
is a significant charge renormalisation effect such that 
Z* < Z by a factor of about 1.2 and 1.5 for Ka = 0.1 and 
Ka = 0.5, respectively. The merging of the red and the 
black curves at high-?/ in FigJIJb) is due to the reduc- 
tion of the charge into the linear- screening regime such 
that Z = Z* . The increase of Z* with observed in 
both FigQa) and (b), is in line with well-known charge- 
renormalisation results, Q22 QH |2"3H21)1 BP] and with Eq. 
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Figure 2: The effective inverse screening length func- 
tion of the packing fraction r\ for several reservoir screening 
parameters na, for constant surface potentials (a) cj>o = 1 and 
(b) (f>o — 5 as represented by the black curves. The red curves 
denote R as obtained from the association-dissociation model, 
with the chargeability z chosen such that the surface poten- 
tial in the dilute limit r\ — > equals 0o- Note that R — k in 
all cases for rj — > 0. 



In Figj2ja) and (b) we plot, for the same zeta- 
potentials as in Fig. 1(a) and (b), the effective screen- 
ing parameter R as a function of 77 for several reservoir 
screening constants n. At low enough 77, where kR 3> 1, 
the two screening constants are indistinguishable from 
each other in all cases. The reason is that the cell is then 
large enough for the potential to decay to essentially zero 
at r = R, such that the asymptotic decay of <fr(r) is gov- 
erned completely by the screening constant k of the back- 
ground (reservoir) salt concentration. At larger 77, and 
hence smaller cells, 4>(R) is no longer vanishingly small 
and the ion concentrations p±(R) at r = R deviate con- 
siderably from the ionic reservoir concentration p s . This 
larger ionic concentration at the cell boundary, which 
represents an enhanced ion concentration in between the 
colloidal particles in the true many-body system, leads 
to a larger effective screening constant k with increasing 
77 at a fixed K, as is shown in Fig{2ja) and (b). Given 
that larger charges are obtained in the AD model than 
in the CSP model at high 77, the number of counterions 
in the cell, and hence R, is also larger in the AD model. 



IV. EFFECTIVE INTERACTIONS AND PHASE 
DIAGRAMS 

Once the effective colloidal charge Z* and the effec- 
tive screening length R~ l have been determined from 
the numerical solution of the PB equation in a Wigner- 
Seitz cell, either for constant-potential or association- 
dissociation boundary conditions, the effective interac- 
tions u(r) between a pair of colloidal particles separated 
by a distance r follows, assuming DLVO theory, as 

, . ( 00, r < 2a; 

^ = I , /Z*exp(Ra)\ 2 exp(-Kr) (12) 

k B T I X B , , _ ' — '-, r>2a\ ' 

\ \ 1 + na J r 

where we include a short-range hard-core repulsion for 
overlapping colloids and ignore Van der Waals forces 
(which is justified for index- matched particles). One 
could use Eq.|l2| for the pair interaction to simulate 
(or otherwise calculate) properties of the suspension in 
a given state-point, e.g. whether the system is in a 
fluid or crystalline state. We restrict our attention 
here to the limiting case in which the colloidal parti- 
cles are sufficiently highly charged and/or sufficiently 
weakly screened, that the pair-potential at contact sat- 
isfies u(2a) ^> ksT, thereby effectively preventing direct 
particle-particle contact. In this limit the suspension can 
be effectively regarded as a point- Yukawa system that 
can be completely characterised by only two dimension- 
less parameters U and A for the strength and the range 
of the interactions, respectively. They are defined as 

= / Z*exp{Ra) \ 2 X B ( 3tA 1/3 

V 1 + Ra ) a \At:) 1 ' 

X = Ra(^y lf \ (14) 

such that the point- Yukawa interaction potential of in- 
terest, in units of ksT, reads [/exp(— Xx)/x with x = 
r(NfV) 1 / 3 the particle separation in units of the typical 
particle spacing. Note that three dimensionless param- 
eters would have been needed if hard-core contact was 
not a low Boltzmann-weight configuration, e.g. then the 
contact-potential f3u(2a) (i.e. the inverse temperature), 
the dimensionless screening parameter na, and the pack- 
ing fraction 77 would be a natural choice. The mapping 
of the phase diagram of the point- Yukawa system onto 
hard-core Yukawa systems has been tested and verified 
explicitly by computer simulation. [14] 



The point-Yukawa system has been studied by simu- 
lation in great detail over the years. [TlTtT4] and by now 
it is well known this model features a disordered fluid 
phase and two crystalline phases (face-centered cubic 
(fee) and body-centered cubic (bcc)). Their first-order 
phase boundaries are well-documented, and can accu- 
rately be described by curves in the two-dimensional 
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Figure 3: Phase diagrams in the packing fraction-screening length (??, « 
a/As = 100) interacting with the hard-core Yukawa potential of Eq.(12l 
lines represent phase boundaries for the constant-potential model, and the red dashed lines for the association-dissociation 
model with the surface potential equal to <f>o in the dilute limit. The dashed black lines indicate extrapolation of Eq. (151 
beyond its strict regime of accuracy. The inset in the phase diagram for (j>o — 5 represents r\ on a logarithmic scale for clarity. 
The labels "Fluid", "BCC", and "FCC" denote the stable fluid, bcc, and fee regions. We note that the very narrow fluid-fee, 
fluid-bee, and fcc-bcc coexistence regions are just represented by single curves. The dotted blue curves represent the estimated 
crossover-packing fraction rf of Eq.|7|, beyond which Z(rj) < Z(0)/2. 



(A, U) plane. Here, we employ the fits for the phase 
boundaries of point-Yukawa particles that were presented 
in Ref . [Hj , which were based on the results of Hamaguchi 
et al. |13j The melting-freezing line between the bcc crys- 
tal and the fluid is accurately fitted by 

\nU = 4.670 -0.04171A + 0.1329A 2 - 0.01043A 3 

+ 4.343 -10" 4 A 4 - 6.924 • 1(T 6 A 5 , (15) 
for < A < 12, 

and the bcc-fcc transition by 

InU = 97.65106 - 150.469699A + 106.626405A 2 
-41.67136A 3 + 9.639931A 4 - 1.3150249A 5 
+0.09784811A 6 - 0.00306396A 7 , (16) 
for 1.85 < A < 6.8. 

Here we exploit these empirical relations as follows. For 
given dimensionless colloid radius a/ Xb and screening 



constant na, we calculate Z* and Ra for various 77 for 
the CSP and the AD model in the Wigner-Seitz cell, as 
described in the previous section. These quantities can be 
used to compute the dimensionless Yukawa parameters U 



and A from Eqs.(13l and (14 1, such that their phase and 



phase-boundaries are known from Eqs.(15l and (16). 



For ci/\b — 100, Figj3] shows the phase diagrams 
that result from this point- Yukawa mapping procedure in 
the (77, (ko)" 1 ) representation, for the CSP model (black 
curves) with surface potentials (a) 0o = 1, (b) <fio = 2, 
(c) 4>o — 3, (d) <j) — 5, and for the corresponding AD 
model (red curves). The dashed lines represent the phase 
boundary fits of Eqs.(15l and (16 1 outside their strict A- 



regime of applicability. We restrict attention to r] < 0.3, 
as the point- Yukawa limit breaks down due to strong 
excluded-volume effects at higher packing fractions. An 
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expected feature is the shift of the freezing curves to 
lower 77 for higher </3o, due to the higher (renormalized) 
charge and the stronger repulsions at higher <f>Q. Due to 
the higher charges in the AD model, its crystallisation 
regimes (red curves) extend to somewhat lower 77's and 
longer screening lengths than those of the CSP model 
(black curves). However, the most striking feature of all 
these phase diagrams is the huge extension of the fluid 
regime: at high and at low screening length there is no 
crystalline phase at all (for 77 < 0.3), while at some inter- 
mediate salt concentrations the crystal phases are sand- 
wiched in between an ordinary low-density fluid and a re- 
entrant fluid phase. This re-entrant fluid regime becomes 
more prominent with increasing zeta-potential 4>q. The 
underlying physics of this finite-salt and finite-77 regime 
where bec and fee crystals exist is the discharging of the 
colloids with increasing 77 and decreasing salt concentra- 
tion: (i) although at high salt (small screening length) 
the colloidal charge is high, the screened-Coulomb inter- 
action is then so short-ranged that the system resem- 
bles a hard-sphere system that will only crystallize at 
77 ~ 0.5; (ii) at low salt (long screening length) the col- 
loidal charge is too low to have sizeable repulsions that 
drive crystallisation. Only at intermediate salt and in- 
termediate colloidal packing the charge is high enough 
and the screening sufficiently long-ranged to drive crys- 
tallisation. The dotted blue curves in Fig [3] represent the 
crossover packing fraction 77* of Eq.Q beyond which the 
colloidal charge has been reduced to less than 50 percent 
of its dilute-limit value. Clearly, our expression for 77* 
indeed roughly coincides with the onset of the re-entrant 
fluid regime. Eq.([7]) thus provides a quick guide to es- 
timate where or whether re-entrant melting is to be ex- 
pected at all. Interestingly, there are parameter values 
for Ka and 4>q (albeit in a narrow range) where a phase 
sequence fluid-bcc-fcc-bcc-fluid is predicted here upon in- 
creasing the colloidal packing fraction, i.e. not only a 
re-entrant fluid phase but also a re-entrant bec phase. 



Moreover, for 77 > 0.5 one expects hard-sphere freezing 
into an fee (or hep) stacking on the basis of hard-sphere 
interactions, so the fee phase is then also re-entrant. 



Experimentally it is not always possible or convenient 
to characterise the screening in terms of the Debye length 
k _1 of the (hypothetical) reservoir with which the sus- 
pension would be in osmotic equilibrium. Rather, one 
often measures the effective (actual) screening length R 
in the suspension of interest. For this reason we replot in 
Fig|4]the phase diagrams for <j> — 2 and <j>o = 5 of Figj3j 
but now in the (77, (Sa) _1 ) representation. Interestingly, 
the CSP and AD model are now much closer together, 
and the re-entrant fluid phase appears even more pro- 
nounced in this representation. 

In order to quantify in which finite salt-concentration 
regime bec and fee crystals are expected in a colloidal 
concentration series < 77 < 0.3, we analyse the maxi- 
mum and minimum values of Ra at which these two crys- 
tal phases can exist, as a function of the zeta-potential 
cj>o, for a/Xs = 100. Fig. 5]shows the resulting screening- 
length regimes, both for bec (black curves) and fee (blue 
curve), where the lowest screening length for fee crys- 
tals is set to zero because of the hard-sphere freezing 
into fee at 77 = 0.5 even for 1/na — > — of course we 
only restricted attention to 77 < 0.3 until now so strictly 
speaking also the fee phase should have had a nonvan- 
ishing lower bound. Nevertheless, despite this small in- 
consistency, Fig(5] clearly shows not only that a larger 
zeta-potential gives rise to a larger crystal regime, but 
also that for all 4>q there is a limiting screening length 
beyond which neither fee nor bec crystals can exist, both 
for the CSP and the AD model. 



So far we focussed on a/Xs — 100, which in aqueous 
suspensions corresponds to a colloidal radius of about 
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Figure 5: Maximum and minimum effective screening lengths where bcc and fee can be found, as a function of the surface 
potential, assuming a constant surface potential for cl/\b = 100. The bcc regime is in between the two black lines, and the fee 
regime below the blue line. The red points indicate the results from the AD model. 
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Figure 6: Phase diagrams in the packing fraction- effective screening length representation (rj, (rea) -1 ), for constant-potential 
colloids with (a) <j>o = 5 for a/Xs = 10 and (b) (j>o = 1 and for q/Xb = 1000. Lines and symbols as in Fig[3] 



70nm. However the colloidal size regime can easily be a 
factor 10 larger or smaller, and for that reason we also 
consider the CSP model for o/A B = 10 and 1000. In Figj6] 
we show the phase diagrams for the smaller colloids with 
</>o = 5 in (a), and for larger colloids with <f) — 1 in (b). 
Interestingly, they resemble those for a/As = 100 shown 
in FigQa) and (b) for </> = 2 and 5, respectively, but the 
smaller colloids require a much higher potential while the 
larger ones only need a lower potential to obtain phase 
diagrams similar to those for a/As — 100 — the phase 
diagram for a/As = 10 at 0o = 1 does not show a crystal 
phase at all for ij < 0.3. In other words, given that 0o = 5 
is a rather high potential that may be difficult to achieve 
in reality while (fto — 1 is frequently occurring, one con- 
cludes that re-entrant melting occurs in the largest salt- 
concentration regime (and is hence easiest observable by 
tuning the salt) for larger colloids. 



V. SUMMARY AND CONCLUSIONS 

Within a Wigner-Seitz cell model we have calculated 
the bare charge Z, the renormalized charge Z* , and the 
effective screening length kT 1 of colloidal spheres at a 
constant zeta-potential cf>o . We find from numerical solu- 
tions of the nonlinear Poisson-Boltzmann equation that 
these constant-potential colloids discharge with increas- 
ing packing fraction and ionic screening length, in fair 
agreement with analytical estimates for the dilute-limit 
charge Zq in Eq.Q and the typical crossover packing 
fraction ?y* given in Eq.([7|. We also show that the 
constant-potential assumption is a reasonably accurate 
description of charge regulation by an ionic association- 
dissociation equilibrium on the colloidal surface. We use 
our nonlinear calculations of Z* and R to determine the 
effective screened-Coulomb interactions between the col- 
loids at a given state point, and we calculate the phase di- 
agram for various zeta-potentials by a mapping onto em- 
pirical fits of simulated phase diagrams of point- Yukawa 



fluids. This reveals a very limited regime of bcc and fee 
crystals: in order to form crystals, the charge is only 
high enough and the repulsions only long-ranged enough 
in a finite intermediate regime of packing fraction and 
salt concentrations; at high r\ or low salt the spheres 
discharge too much, and at high salt the repulsions are 
too short-ranged to stabilise crystals. In the salt-regime 
where crystals can exist, the discharging mechanism gives 
rise to re-entrant phase behaviour, with phase sequences 
fluid-bec-fluid and even fluid-bcc-fcc-bcc-fluid upon in- 
creasing the colloid concentration from extremely dilute 
to 77 = 0.3. 

The phase behaviour of constant-potential or charge- 
regulated colloids as reported here is quite different from 
that of constant-charge colloids, for which the pairwise 
repulsions do not weaken with increasing volume frac- 
tion or decreasing salt concentration. As a consequence 
constant-charge colloids have a much larger parameter- 
regime where crystals exist, and do not show the re- 
entrant behaviour. |lllfl"4] The most direct comparison 
is to be made with the constant-charge phase diagrams 
of Fig. 2 and Fig. 4 of Ref. [14] . where the charge is fixed 
such that the surface potential at infinite dilution cor- 
responds to 4>o — 1 and 2, respectively. Our theoretical 
findings can thus be used to gain insight into the colloidal 
charging mechanism by studying colloidal crystallisation 
regimes as a function of packing fraction and salt con- 
centration. 
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Appendix 

Although it is numerically straightforward to solve the 
nonlinear PB equation ^ with BC's ^ and Q in a 
spherical Wigner-Seitz cell of radius R, it may also be 
convenient to have analytic results that allow for quick es- 
timates of the (order of) magnitude of the colloidal charge 
Z. A standard approach is to linearise the smb(f>(r) term 
of Eq.|2]), e.g. with <j)(r) — 4>{R) as the small expansion 
parameter. The resulting solution is then of the form 
4>(r) — Aenp(—Rr)/r + B exp(Rr) / V + C, with R defined 
in Eq.(ll|, C = <fi(R) — tanh(/>(i?), and with integration 
constants A and B fixed by the two BC's. The algebra 
involved is, however, not very transparent. 

A considerable simplification is achieved if we consider 
the so-called Jellium model, in which the central col- 



loidal sphere is no longer considered to be surrounded 
by only cations and anions in a finite cell, but instead by 
cations, anions and other colloids with charge Z (to be 
determined) .[24 26 A nonlinear PB equation and BC's 
can then be written, for r > a. 



0"(r) + -0'(r) 
r 

<t>{a) 
<P'(oo) 



k 2 smb. <j)(r) - AttAbZu; (17) 



0. 



(18) 
(19) 



where it is assumed that the 'other' colloids are dis- 
tributed homogeneously with density n. From this one 
derives directly that the asymptotic potential is given by 



AttXeZh 3r/(ZX B /a) 
sinh</>(oo) = = 



(20) 



Now linearising sinh 4>{r) with 4>(r) — </>(oo) as the small 
expansion parameter gives rise to the solution 

cf>(r) = 0(oo) + (0o - 0(oo)) 6XP ( - K(r - a)) , (21) 

r/a 

where the effective screening length k^ 1 is defined by 

V^osh^oo). (22) 

We note that the average ion concentrations in the sys- 
tem, within the present linearisation scheme, is given 
by c± = p s exp(=F0(oo)), such that the corresponding 
screening length k^ 1 is given by k 2 = / LttXb(c+ + c_). 
In other words, the effective screening length k and the 
asymptotic potential 0(oo) of this jellium model play ex- 
actly the same role as R and <fr(R) that we introduced 
before in the Wigner-Seitz cell. In particular, R^ 1 and 
k _1 can be seen as the actual screening length in the sus- 
pension (in contrast to the screening length k _1 of the 
ion reservoir). 



From Eq.(21) the colloidal charge Z follows, using 
Eq. |5]) , as the solution of the transcendental equation 



ZX f 



= (0o ~ 0(oo))(l + ka), 



(23) 



where one should realise that both (j)(oo) and k depend 
on ZXs/a through Eqs.(20) and (22). It is possible to 
solve Eq.(23) explicitly in the dilute limit. For 77 = 
one finds 



00 



from Eq.((20f, and hence Z = Z a 
given by Eq.Q. For finite but low-enough rj for which 
</>(oo) <C 1 one can ignore 0(r/ 2 ) contributions, such that 
sinh (f>(oo) ~ 0(oo) and cosh0(oo) ~ 1, to find Eq.([7]) 



from the self-consistency condition Eq. ( 23 1 
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